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Triplet superconductivity, which is an analogue of su- 
perfluidity in 3 He, was firstly discovered in a heavy- 
fermion compound UPt^ more than ten years ago 
Recently, a novel triplet superconductor 5Y2-RUO4 was 
found In these compounds, the triplet pairing can 

be reliably determined, for example, by Knight shift ex- 
periments , however the identification of a symmetry 
of the order parameter is much more difficult task. A 
large number of experimental and theoretical investiga- 
tions done on UPt^ and S^RuOi are concerned with 
different thermodynamic and transport properties, but 
the precise order parameter symmetry remains still to be 
worked out yet (see, for example, [t|J10|-[i^| , and original 
references therein). 

Calculations of the order parameter A in U Pts 
and Sr2RuO^ as a function of the momentum direction 
k on the Fermi surface is a very complex problem. Some 
general information about A can be obtained from 
a symmetry of a normal state: G sp i n - or i>it x r x U(l), 
where G SV i n -orbit represents the point group with inver- 



A stationary Josephson effect in point contacts be- 
tween triplet superconductors is analyzed theoretically for 
most probable models of the order parameter in VPtz and 
Sr2RuOi. The consequence of misorientation of crystals in 
superconducting banks on this effect is considered. We show 
that different models for the order parameter lead to quite dif- 
ferent current-phase dependences. For certain angles of mis- 
orientation a boundary between superconductors can generate 
the parallel to surface spontaneous current. In a number of 
cases the state with a zero Josephson current and minimum of 
the free energy corresponds to a spontaneous phase difference. 
This phase difference depends on the misorientation angle and 
may possess any value. We conclude that experimental inves- 
tigations of the current-phase dependences of small junctions 
can be used for determination of the order parameter symme- 
try in the mentioned above superconductors. 

I. INTRODUCTION. 



sion; t is the time-inversion operator, and U(l) is the 
group of gauge transformation. A superconducting state 
breaks one or more symmetries. In particular, a tran- 
sition to the superconducting state implies an appear- 
ance of the phase coherence corresponding to breaking of 
the gauge symmetry. According to Landau theory fjlf l 
of second order phase transitions, the order parameter 
is transformed only on irreducible representations of the 
symmetry group of the normal state. Conventional su- 
perconducting states have a total point symmetry of the 
crystal and belong to the unitary even representation 
Ai g . In unconventional superconductors this symmetry 
is broken. The parity of a superconductor with inver- 
sion symmetry can be specified using the Pauli principle. 
Because for triplet pairing the spin part of the A is a sym- 
metric second rank spinor, the orbital part has to belong 
to the odd representation. In a general case the triplet 
paring is described by the order parameter of the form 
A ^k^ = id er<72, where the vector a — (i^i, 02, (73) , 
and Ui are Pauli matrices in the spin space. A vector 
d = — d ^— k^j in spin space is frequently referred 
to as an order parameter or a gap vector of the triplet 
superconductor. This vector defines the axis along which 
the Cooper pairs have zero spin projection. If d is com- 
plex, the spin components of the order parameter spon- 
taneously break time-reversal symmetry. 

Symmetry considerations reserves for the order param- 
eter considerable freedom in the selection of irreducible 
representation and its basis functions. Therefore in many 



papers (see, for example, [[7pt^- |T^Jl^ |l6| |) authors con- 
sider different models (so-called scenarios) of supercon- 
ductivity in UPtz and Sr2RuOi, which are based on 
possible representations of crystallographic point groups. 
The subsequent comparison of theoretical results with 
experimental data makes it possible to conclude on the 
symmetry of the order parameter. 

In real crystalline superconductors there is no classifi- 
cation of Cooper pairing by angular momentum (s-wave, 
p-wave, d-wave, /-wave pairing, etc.). However these 
terms are often used for unconventional superconductors 
meaning that the point symmetry of the order parameter 
is the same as one for the corresponding representation of 
SO3 symmetry group of isotropic conductor. In this ter- 
minology conventional superconductors can be referred 
to as s-wave. For example, the "p-wave" pairing corre- 
sponds to the odd two-dimensional representation E\ u of 
Deh point group or E u representation olD^h point group. 
The order parameter for these representations has the 
same symmetry, as for the superconducting state with 
angular momentum I = 1 of Cooper pairs in an isotropic 
conductor. If the symmetry of A can not be formally 
related to any irreducible representation of SO3 group, 
these states are usually referred to as hybrid states. 

Apparently, in crystalline triplet superconductors the 
order parameter has more complex dependence on k in 
comparison with well known p-wave order parameter for 
superfluid phases of 3 He. The heavy-fermion supercon- 
ductor UPtz belongs to the hexagonal crystallographic 



point group (Dq^) > and it is most likely that the pairing 
state belongs to Ei u ("/—wave" state) representation. 
A layered perovskite material Sr^RuO^ belongs to the 
tetragonal crystallographic point group (D^) . Initially 
the simplest "p-wave" model based on E u representation 
was proposed for the superconducting state in this com- 
pound 1^,^]. However this model was inconsistent with 
available experimental data, and later [|lO],[ll| other "/- 
wave" models of pairing state were proposed. 

Theoretical studies of specific heat, thermal conductiv- 
ity, ultrasound absorption for different models of triplet 
superconductivity show considerable quantitative differ- 
ences between calculated dependences Q|i^,[ll],|l6| . The 
Josephson effect is much more sensitive to dependence 
of A on the momentum direction on the Fermi surface. 
One of the possibilities to form a Josephson junction is to 
create a point contact between two massive superconduc- 
tors. A microscopic theory of the stationary Josephson 
effect in ballistic point contacts between conventional su- 
perconductors was developed in Ref. |l^| . Later this the- 
ory was generalized for a pinhole model in 3 He [ ^8|Jl9| 
and for point contacts between "d-wave" high-T c super- 
conductors |2^j2l[]. It was shown that current-phase de- 
pendences for the Josephson current in such systems are 
quite different from those of conventional superconduc- 
tors, and states with a spontaneous phase difference be- 
come possible. Theoretical and experimental investiga- 
tions of this effect in novel triplet superconductors seem 
to be interesting and enable one to distinguish among 



different candidates for the superconducting state. 

In Ref. |2^| the authors study the interface Andreev 
bound states and their influence on the Josephson cur- 
rent between clean "/-wave" superconductors both self- 
consistently (numerically) and non-self-consistently (an- 
alytically). The temperature dependence of the critical 
current is presented. However in that paper there is no 
detailed analysis of the current-phase dependences for 
different orientations of the crystals in the superconduct- 
ing banks. 

In this paper we theoretically investigate the stationary 
Josephson effect in a small ballistic junction between two 
bulk triplet superconductors with different orientations 
of the crystallographic axes with respect to the junction 
normal. In Sec. II we describe our model of the junction 
and present the full set of equations. In Sec. Ill the cur- 
rent density in the junction plane is analytically calcu- 
lated for a non-self-consistent model of the order param- 
eter. In Sec. IV the current-phase dependences for most 
likely models of " /-wave" superconductivity in U Pt^ and 
S^RuO^ are analyzed for different mutual orientations 
of the banks. We end in Sec.V with some conclusions. 

II. MODEL OF THE CONTACT AND 
FORMULATION OF THE PROBLEM. 

We consider a model of a ballistic point contact as 
an orifice with a diameter d in an impenetrable for elec- 
trons partition between two superconducting half spaces 
(Fig.I). We assume that the contact diameter d is much 
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FIG. 1. Scheme of a contact in the form of an orifice be- 
tween two superconducting banks, which are misorientated 
on an angle a. 

larger than the Fermi wavelength and use the quasiclassi- 
cal approach. In order to calculate the stationary Joseph- 
son current in point contact we use " transport-like" equa- 
tions for ^-integrated Green functions 9 ^k, r, 23 



i£ m t 3 - A, 9 



iv F kV .9= 0, 



and the normalization condition 



vv 
.9.9= 



-1. 



(1) 



(2) 



Here e m — nT(2m + 1) are discrete Matsubara energies, 
vp is the Fermi velocity, k is a unit vector along the elec- 
tron velocity, t 3 — T3<g> J; (i = 1, 2, 3) are Pauli matrices 
in a particle-hole space. 

v 

The Matsubara propagator 9 can be written in the 
form g§]: 

/ 



9i + gi^ (.92 + g2<?) ia 2 

\ i&2 (.93 + g3<?) 9i - V2g4Ve2 



(3) 



as can be done for an arbitrary Nambu matrix. Ma- 

v 

trix structure of the off-diagonal self energy A in Nambu 
space is 



v 



(4) 



id(7(72 

l id 2 d*d 

Below we consider so-called unitary states, for which 
d x d* = 0. 

The gap vector d has to be determined from the self- 
consistency equation: 

d (k, r) = ttTN (0) (V (k, k') g 2 (k', r, e m ) ) , (5) 

rn 

where V ^k, k'^ is a potential of pairing interaction; (...) 
stands for averaging over directions of an electron mo- 
mentum on the Fermi surface; N (0) is the electron den- 
sity of states. 

Solutions of Eqs. (@), (§ must satisfy the conditions 
for Green functions and vector d in the banks of super- 
conductors far from the orifice: 



v v 

X , N Km ^3 - Al,2 

9 (Too) = 



|di. 



d(=Foo) = di 2 (k) exp ( =f 



icj) 



(6) 
(7) 

where <fi is the external phase difference. Eqs. (|l|) and 
(||) have to be supplemented by the boundary continuity 
conditions at the contact plane and conditions of reflec- 
tion at the interface between superconductors. Below we 
assume that this interface is smooth and electron scat- 
tering is negligible. 

III. CALCULATION OF THE CURRENT 
DENSITY. 

The solution of Eqs. (|l|) and (||) allows us to calculate 
the current density: 



j (r) = 2neTv F N (0) ]T (kgi (k, r, e m ) ) 



(8) 



We consider a simple model of the constant order pa- 
rameter up to the surface. The pair breaking and the 
scattering on the partition and in the junction are ig- 
nored. This model can be rigorously found for the calcu- 
lations of the current density (||) in ballistic point con- 
tacts between conventional superconductors at zero ap- 
proximation with small parameter d/£o (£o is the coher- 
ence length) jnj . In anisotropically paired superconduc- 
tors the order parameter changes at distances of the or- 
der of £o even near a specular surface Thus for 
calculations of the current (^|) in the main approxima- 
tion on the parameter d/£o it is necessary to solve the 
equation (||) near a surface of the semi-infinite super- 
conductor. It can be done only numerically and will be 
the subject of our future investigations. Below we as- 
sume that the order parameter does not depend on co- 
ordinates and in each half-space equals to its value (0) 
far form the point contact. For this non-self-consistent 
model the current-phase dependence of a Joscphson junc- 
tion can be calculated analytically. It makes possible to 
analyze the main features of current-phase dependences 
for different scenarios of " /-wave" superconductivity. We 
believe that under this strong assumption our results de- 
scribe the real situation qualitatively, as it was justified 
for point contacts between " d- wave" superconductors 
and pinholes in 3 He [^tJ. It was also shown in Ref. |2|] 
that for the contact between "/-wave" superconductors 
there is also good qualitative agreement between self- 



consistent and non-self-consistent solutions (although, of 
course, quantitative distinctions, take place). 

In a ballistic case the system of 16 equations for func- 
tions gi and gi can be decomposed on independent blocks 
of equations. The set of equations which enables us to 
find the Green function g\ is 

iv F kVgi + (g 3 d - g 2 d*) = 0; (9) 

w F kVg_ + 2i (d x g 3 + d* xg 2 ) = 0; (10) 

iu F kVg 3 - 2ie m g 3 - 2gid* - id* x g_ = 0; (11) 

iv F kVg 2 + 2ie m g 2 + 2gid - id x g_ = 0; (12) 

where g- = gi — g4- The Eqs. (|9|)-(|l2}) can be solved by 
integrating over ballistic trajectories of electrons in the 
right and left half-spaces. The general solution satisfying 
the boundary conditions at infinity is 



g W = ^ + lCn exp (-2sQ n t) ; (13) 
= C n exp (-2sfi„t) ; (14) 

^"- 2 -t^^" ° p( - 2 ' n "' ) -^ ; (15) 

g („, = 2c„d; + d;x c„ cxp ( _ 2<fu) _ d; . (M) 

where t is the time of the flight along the trajectory, 



sgn (t) = sgn (z) = s; r\ = sgn (v z ) ; ft n = \j e 2 m + |d„| 2 . 
By matching the solutions (p^-^6[) at the orifice plane 
(t = 0), we find constants C n and C n . Index n numbers 
left (n = 1) and right (n = 2) half-spaces. The function 
gi (0) = g^ (—0) = <?i (+0) , which determines the cur- 
rent density in the contact is 

9i (0) = 



ie m cose (^i + fi 2 ) + fysine (e^ + f2iSl 2 ) 



(17) 



AiA 2 + (e 2 „ + f2i0 2 ) cose - is m V (^1 + fi 2 ) sine 
In the formula ( |l7| ) we have taken into account that for 
unitary states vectors di 2 can be written as 



d„ = A„exp itp n , 



(18) 



where Ai 2 are real vectors. 

Knowing the function g\ (0) one can calculate the cur- 
rent density at the orifice plane j (0) : 

00 

j(0)=4ireN{0)v F TY^ dkkRe gi (0), (19) 



m=0 



where 



Re 9l (0) = 

sine (A|A| cose + (e^ 4- »i^ 2 ) AiA 2 ) 
[Ax A 2 + ( e 2, + fi^a) cose] 2 + (fli + ft 2 ) 2 sin 2 e 



(20) 



or, alternatively, 



Re ffl (0) 



AiA 2 x ^ 



sin(e ± ' 



2 <^ e 2 ,, + + A^cos^i!?)' 

(21) 



where is defined by Ai(k)A 2 (k) = A x (k)A 2 (k) cos6>, 
ande(k)=V2(k)-Vi(k) + 0. 

Misorientation of the crystals would generally result 
in appearance of the current along the interface [ ^0|j22| 
as can be calculated by projecting vector j at the corre- 
sponding direction. 

We consider a rotation R only in the right supercon- 
ductor (see, Fig.l), (i.e., d 2 (k) = itai (i^k)). The 
c-axis in the left half-space we choose along the partition 
between superconductors (along z-axis in Fig.l). To il- 



we plot the current-phase dependence for different below 
mentioned scenarios of " /-wave" superconductivity for 
two different geometries corresponding to different orien- 
tations of the crystals to the right and to the left at the 
interface (see, Fig.l): 

(i) The basal afr-plane to the right is rotated about 
c-axis by the angle a; Ci||c 2 . 

(ii) The c-axis to the right is rotated about the contact 
axis (y-axis in Fig.l) by the angle a; bi||b 2 . 

Further calculations require a certain model of the vec- 
tor order parameter d. 

IV. CURRENT-PHASE DEPENDENCE FOR 
DIFFERENT SCENARIOS OF "F-WAVE" 
SUPERCONDUCTIVITY. 

The model which has been successful to explain prop- 
erties of the superconducting phases in UPt^ is based 
on the odd-parity Ei u representation of the hexago- 
nal point group for the strong spin-orbital cou- 
pling with vector d locked along c axis of the lattice 
@: d =A z[r ll Y 1 + m Y 2 ] , where Y x = k z (fc 2 - fc 2 ) and 
Y2 = 2fc x fcj / fc z are the basis function of the representa- 
tionPJ The coordinate axes x,y,z here and below are 
chosen along the crystallographic axes a, b,c as at the 



lustrate results obtained by computing the formula (19) 



1 Strictly speaking, in crystals with a strong spin-orbit cou- 
pling the spin is a "bad" quantum number, but electron 
states are twice degenerated and can be characterized by 
pseudospins. 



left at Fig.l. This model describes the hexagonal analog 
of spin-triplet /-wave pairing. For the high-temperature 
A-phase (772 = 0) the order parameter has an equato- 
rial line node and two longitudinal line nodes. In the 
low-temperature B-phase (772 = i) or the axial state: 

d =A zk z (k x + iky) 2 , (22) 

the longitudinal line nodes are closed and there is a 
pair of point nodes. The B-phasc (|22] ) breaks the time- 
reversal symmetry. The function Ao = Ao (T) in Eq.|22] 
and below describes the dependence of the order param- 
eter d on temperature T (carrying out numerical calcu- 
lations we assume T = 0). 

Other candidates to describe the orbital states, which 
imply that the effective spin-orbital coupling in UPtj, is 
weak, are unitary planar state 

d=A k z (x (kl-k 2 .) +y2k x k y ), (23) 

(or d =Ao (Yi, I2, 0)) and non-unitary bipolar state 
d =A (Yi,iY 2 ,0) 0. In Fig.| we plot the Josephson 
current-phase dependence jj(4>) — j y (y — 0) calculated 
from Eqjl^ for both the axial (with the order parameter 
given by Eq.^2|) and the planar (Eqj2^) states for a par- 
ticular value of a under the rotation of basal afr-plane to 
the right (the geometry (i)). For simplicity we use the 
spherical model of the Fermi surface. For the axial state 
the current-phase dependence is just the slanted sinusoid 
and for the planar state it shows a "7r-state". The ap- 
pearance of 7r-state at low temperatures is due to the fact 
that different quasiparticle trajectories contribute to the 



■ planar 




FIG. 2. Josephson current densities versus phase <f> for ax- 
ial (22) and planar (23) states in the geometry (i); misori- 
entation angle a — n/4; the current is given in units of 
jo = f eN(0)v F A o (0). 

current with different effective phase differences ?(k) (see 
Eqs.|l9| and |2l]) [jl9). Such a different behavior can be a 
criterion to distinguish between the axial and the planar 
states, taking advantage of the phase-sensitive Joseph- 
son effect. Note that for the axial model the Josephson 
current formally does not equal to zero at <f> = 0. This 
state is unstable (does not correspond to a minimum of 
the Josephson energy), and the state with a spontaneous 
phase difference (value </>o in Fig. 2), which depends on 
the misorientation angle a, is realized. 

The remarkable influence of the misorientation angle a 
on the current-phase dependence is shown in FigJ^ for the 
axial state in the geometry (ii) . For some values of a (in 
Fig|| it is a = 7r/3 ) there are more than one state, which 
correspond to minima of the Josephson energy (jj = 
and djj/dxj) > 0). 

Calculated x and z-components of the current, which 



FIG. 3. Josephson current versus phase <jj f° r the axial (22) 
state in the geometry (ii) for different a. 

are parallel to the surface, js(</>) are shown in Fig.[| and 
Fig.|| for the same axial state in the geometry (ii) . Note 
that the tangential to the surface current as a function 
of (p is not zero when the Josephson current (Fig||) is 
zero. This spontaneous tangential current (see also in 
Rcf. |2^|) is due to the specific "proximity effect" simi- 
lar to spontaneous current in contacts between " c£-wave" 



superconductors 1 20 2q] . The total current is determined 
by the Green function, which depends on the order pa- 
rameters in both superconductors. As a result of this, 
for nonzero misorientation angles the parallel to the sur- 
face current can be generated. In the geometry (i) the 
tangential current for both the axial and planar states at 
T = is absent. 

The first candidate for the superconducting state in 
S^RuOi was "p-wave" model [p| 



z (k x + iky^j . 



d =A„ : f /■■., + //■■„) . (24) 

Recently Jll|,|12| it was shown that the pairing state in 
5r2-RwC?4, most likely, has lines of nodes. It was sug- 
gested that this can occur if the spin-triplet state be- 



FIG. 4. Z-component of the tangential current versus phase 
for the axial state (22)in the geometry (ii) for different a. 



FIG. 5. X-component of the tangential current versus 
phase (j) for the axial state (22)in the geometry (ii) for dif- 
ferent a. 

longs to a non trivial realization of the E u representation 
of D 4h group, either B 2g <8> E u @ or B lg ® E u @ sym- 
metry: 

d =Ao'zk x k y (k x + ik y ^j , for B^g ® E u symmetry; 

(25) 



d =A z [k x — ky) ( k x + iky ) , for B\ g ® E u symmetry. 



(26) 



Note that models (E4-Era)of the order parameter sponta- 



FIG. 6. Josephson current versus phase <f> for hybrid 
"/-wave" and "p-wave" states in the geometry (i); a = 7r/4. 

neously break time-reversal symmetry. 

Taking into account a quasi-two-dimensional electron 
energy spectrum in Sr2RuOi, we calculate the current 
( |l9| ) numerically using a model of the cylindrical Fermi 
surface. The Josephson current for the hybrid "/-wave" 
model of the order parameter Eq.(p6|) is compared to 
the p-wave model (Eq.^4|) in Fig (for a — 7r/4). Note 
that the critical current for the " /-wave" model is several 
times smaller (for the same value of Ao) than for the 
"p-wave" model. This different character of the current- 
phase dependencies enables us to distinguish between the 
two states. 

In Figs. and|| we present the Josephson current and 
the tangential current for the hybrid " /-wave" model for 
different misorientation angles a (for the "p-wave" model 
it equals to the zero). Just as in Fig. 2 for the "/-wave" 
order parameter (22), in Fig. 7 for the hybrid "/-wave" 
model (25) the steady state of the junction with zero 
Josephson current corresponds to the nonzero sponta- 



FIG. 7. Josephson current versus phase <f> for the hybrid 
"/-wave" state in the geometry (i) for different a. 
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FIG. 8. Tangential current density versus phase <f> for the 
hybrid "/-wave" state in the geometry (i) for different a. 

neous phase difference if misorientation angle a =/= 0. 

V. CONCLUSION. 

We have considered the stationary Josephson effect 
in point-contacts between triplet superconductors. Our 
consideration is based on models with "/-wave" sym- 
metry of the order parameter belonging to the two- 
dimensional representations of the crystallographic sym- 
metry groups. It is shown that the current-phase de- 
pendences are quite different for different models of the 



order parameter. Because the order parameter phase de- 
pends on the momentum direction on the Fermi surface, 
the misorientation of the superconductors leads to spon- 
taneous phase difference that corresponds to the zero 
Josephson current and to the minimum of the weak link 
energy. This phase difference depends on the misorien- 
tation angle and can possess any values. We have found 
that in contrast to "p-wave" model, in the "/-wave" mod- 
els the spontaneous current may be generated in a direc- 
tion which is tangential to the orifice plane. Generally 
speaking this current is not equal to zero in the absence 
of the Josephson current. We demonstrate that the study 
of a current-phase dependence of small Josephson junc- 
tion for different crystallographic orientations of banks 
enables one to judge on the applicability of different mod- 
els to the triplet superconductors UPt 3 and Sr 2 Ru04. 

It is clear that such experiments require very clean su- 
perconductors and perfect structures of the junction be- 
cause of pairbreaking effects of non-magnetic impurities 
and defects in triplet superconductors. The influence of 
single impurities and roughness of interface in the plane 
of the contact, which may essentially decrease the criti- 
cal current of the junction, will be analyzed in our next 
paper. 
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